Fortschritte der Physik, 4 March 2013 



Neutrino Mixings and the 5 4 Discrete Flavour Symmetry 

Federica Bazzocchi 1 and Luca Merlo 2 3 * 

1 INFN - Sezione di Trieste and SISSA, via Bonomea 265, 34136 Trieste, Italy, 

2 Physik-Department, Technische Universitat Miinchen, James-Franck-Strasse, D-85748 Garching, Ger- 
many-preprint: TUM-HEP-834/12 

3 TUM Institute for Advanced Study, Technische Universitat Miinchen, Lichtenbergstrasse 2a, D-85748 
Garching, Germany 

Received XXXX, revised XXXX, accepted XXXX 
Published online XXXX 

Key words Discrete Symmetries, Neutrino Mixings, Lepton Flavour Violation 

Discrete non-Abelian Symmetries have been extensively used to reproduce the lepton mixings. In particu- 
lar, the S4 group turned out to be suitable to describe predictive mixing patterns, such as the well-known 
Tri-Bimaximal and the Bimaximal schemes, which all represent possible first approximations of the ex- 
perimental lepton mixing matrix. We review the main application of the S4 discrete group as a flavour 
symmetry, first dealing with the formalism and later with the phenomenological implications. In particular, 
we summarize the main features of flavour models based on S4, commenting on their ability in reproduc- 
ing a reactor angle in agreement with the recent data and on their predictions for lepton flavour violating 
transitions. 
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1 Introduction 

The presence of additional symmetries beyond the SM gauge group acting on the three fermion generations 
can produce realistic mass hierarchies and mixing textures. In such models, the Lagrangian is invariant 
under the gauge group of the SM and under the additional flavour symmetry at an energy scale A f equal 
or higher than the electroweak one. Fermion masses and mixings then arise once these symmetries are 
broken, spontaneously or explicitly. These flavour models differ from each other depending on the nature 
of the symmetries, either global or local, either Abelian or non-Abelian, either continuous or discrete, and 
on their symmetry breaking mechanism. Many examples have been produced in the last thirty years: at 
first only the quark sector has been investigated and by the implementation of the simple U(l) ITj and 
U(2) (2}|4j symmetries; more recently, also the lepton sector received lot of attention, due to the more 
precise data on the neutrino oscillations |5]|6). In particular, discrete non-Abelian symmetries |[7]-[T2| have 
shown their ability in describing particular mixing patterns, that well describe the experimental data: these 
textures are the so-called Tri-Bimaximal p3||T7l (TB), Bimaximal (T8}^TJ (BM), Golden Ratio (22}|27) 
(GR) and Trimaxinal |28|(TM) schemes. 

Apart from the symmetries, flavour models differ from each other depending on the choice of the flavour 
symmetry breaking scale, that could either correspond to the electroweak scale or be larger, up to the 
grand unification scale. In the first case, replicants of the SM Higgs, with non-trivially transformation 
properties under the flavour symmetry, usually enrich the spectrum. These models could suffer of strong 
constraints from Flavour Changing Neutral Current (FCNC) processes, that receive contributions from 
Higgs exchanges [29-34]. In the second case, new scalar fields, invariant under the SM gauge group but 



transforming under the flavour symmetry, enrich the spectrum of the models. Their masses and vacuum 
expectation values (VEVs) are typically very large, close to the grand unification energy scale, and as 
a result they have negligible direct contributions to FCNC processes. Indirect contributions, however, 
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are possible in beyond SM contexts with the presence of new physics (NP) |35|, such as supersymmetry 
|36H46) or extra dimensions (ED) (47]|48). 

In order to avoid large contributions to FCNC processes, it could be useful to follow a bottom-up 
approach: it consists first in identifying low-energy effective schemes in which the contributions to FCNC 
observables are under control and subsequently in constructing high-energy models which could origin 
such effective descriptions. The so-called Minimal Flavour Violation scheme (MFV) [49 50] follows this 
approach. Since no evident (larger than 3er's) deviations from the SM predictions have been found in all 
flavour processes observed in the hadronic sector, from rare decays in the kaon and pion sectors to B 
decays at superB-factories fST) , it seems reasonable that any physics beyond the SM does not introduce 
new sources of flavour and CP violation with respect to the SM. In Refs. |52-56| this criterion has been 
rigorously defined in terms of continuous flavour symmetries, considering an effective operator description 
in the SM context. Even if the MFV ansatz is successful to prevent large flavour violation also in models 
beyond SM J57fj63] l, it does not provide any explanation for the origin of fermion masses and mixings. 
From this the idea that the MFV is a too restrictive context to provide any clue on the origin of masses. 
Some extensions of the MFV where smaller continuous symmetries have been adopted can be found in 
Refs. J64H661, but also in these analyses a complete explanation of fermion masses is missing. 

In this review we deal with discrete non-Abelian symmetries that lead to a more complete explanation 
of the fermion masses. Furthermore, we restrict our analysis only on models where the flavour symmetry 
is broken at a high scale, that is a more promising framework to prevent large flavour violation. More 
in details, we concentrate on the S4 discrete group that has been adopted to describe the TB, BM and 
TM neutrino patterns. In Sect. [2] we motivate the choice of discrete symmetries to handle the neutrino 
sector and describe the aforementioned mixing schemes, commenting on their agreement with the data. In 
Sect. [3] we describe the ^4 group and three relevant bases, that have been used in the literature. In Sect. |4| 
we summarize the main features of flavour models dealing with the S4 group, including the extension to 
the quark sector. In Sect. [5] we study some phenomenological implications of the S4 models, such as 
the presence of neutrino mass sum-rules, the prediction for the neutrino-less-double-beta (0z/2/3) decay 
effective mass and the constraints coming from lepton flavour violating decays. Finally, in Sect. [6] we 
conclude. 



2 Lepton Mixing Patterns 

The precision data era for lepton mixing started only ten years ago but since then our knowledge of lepton 
mixing angles as well as that of neutrino mass splitting have grown exponentially. Nowadays, the most 
recent results on the oscillation data analysis^can be summarized in Tab. [I] 



Table 1 Recent fits to neutrino oscillation data from (5]|6J 
full (0, 2tt) is allowed. 
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Notice that the best fit values for the reactor angle differ for a factor of 2 in the two fits, due to the exclusion (Fogli et al.) or 
the inclusion (Schwetz et al.) of the data from SBL neutrino experiments with a baseline < 100 m. 
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The past and present years were important for neutrino physics because in the summer 2011 the last 
T2K data [67 1 showed evidences for a non-vanishing reactor angle at the 3<7 level. Subsequently also 
MINOS [68 1 and Double Chooz J69) presented their results, in agreement with T2K one. More recently the 
Daya Bay [70] and RENO fTT) experiments have released their results on the observation of electron anti- 
neutrino disappearance, providing at more than 5<7 and 6a, respectively, the evidence for a non-vanishing 
reactor angle: 

sin 2 6»i3 = 0.024 ±0.005 [Daya Bay], sin 2 13 = 0.029 ±0.006 [RENO] . (1) 

The average of the three most precise results for the reactor angle from DOUBLE CHOOZ (69), Daya 
Bay ff0\ and RENO [71] experiments is given by 

sin 2 0i3 = 0.0253 ± 0.0035 , (2) 



for both the mass orderings. 

These results indicate the double nature of neutrino puzzle: indeed, the neutrino mass scale and mass 
hierarchies are smaller of many orders of magnitude with respect to those of the other fermions and the 
lepton mixing is completely different from quark mixing, where the largest mixing angle is the Cabibbo 
angle A = sin 6^ — 0.22. For these reasons neutrino physics has raised so much interest in the last 
years. In this review we will focus on the origin of lepton mixing and will not address the other interesting 
question related to the origin of their masses. In all the models we will present neutrino masses are obtained 
by means of See-Saw mechanisms or more simply by means of the Weinberg operator. 

Without any kind of prejudice we may state that there are two motivations in support of the idea that 
lepton mixing is ascribed to the breaking of non-Abelian discrete flavour symmetry. Firstly, the presence 
of two large angles have suggested the idea that the lepton mixing matrix cannot be well described by mass 
ratios, as for the CKM matrix. In order to encode this ansatz, in Ref. f72| the concept of Form Diagonal- 
izable (FD) matrices has been introduced: that is mass matrices that are diagonalized by unitary matrices 
whose entries are independent of the mass eigenvalues. In the following, with a small abuse of notation, we 
will refer to these unitary matrices as FD mixing matrices even if the FD acronym refers only to the mass 
matrices. Secondly, under the assumption that neutrinos are majorana particles, if there is any residual 
symmetry behind neutrino mass matrix this is at most a Z2 x Z 2 flavour symmetry j 73 74). The technical 
implication of this feature will be seen in next sections. In the rest of the section, we concentrate on the 
predictive FD mixing patterns. 



First of all we remind the standard parametrization adopted for lepton mixing. This is written in terms 
of 3 mixing angles, 1 Dirac phases and 2 Majorana phases. The three mixing angles, 612, 023, 613 are also 
denoted as solar, atmospheric and reactor angles in order to remember the kind of neutrino source. The 
standard form is given by 

Ul ep = i?23(^23) ' Ria(9l3,S) ' -^12(^12) ' P(f> 




10 \ / cis s l3 e iS \ ( c 12 an \ (3) 

C 23 S 23 • 1 • S U C 12 • P 

S23 C23 J \ -Sl3e~ lS C13 / V 1 / 



with RijiOij) a rotation in the i — j plane of an angle dij,P§ = diag(l, e"* 21 , e M31 ), where ctk are the Ma- 
jorana phases, and c,j , Sij = cos 6ij , sin 6ij . According to Ref. J72), a FD mixing pattern is characterized 
by having all the Cy, as pure numbers, in particular independent from the mass eigenstates. For this 
reason a specific FD mixing pattern is not correlated neither to the neutrino spectrum nor to the Majorana 
phases. In concrete models the constraints on the neutrino mass splittings typically provide specific ranges 
for the Majorana phases. More these ranges are small more predictive, and thus interesting, are the models. 
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The first FD mixing pattern that we consider was suggested by Harrison Perkins Scott in 2002 1 13 T7J 
and is now referred as the Tri-Bimaximal (TB) scheme: 



Utb 



( -2- - 1 - 

1_ J_ 1_ 

\ V6 \/3 V2 



(4) 



that corresponds to the following mixing angles 

sin 023 = \ , sin Of 2 = ^ , sin0? 3 = . (5) 
The most general mass matrix diagonalized by the TB mixing is given by 

x y y 
y x + z y — 
y y — z x+z 



m l B = \ V x + z y-z \ , (6) 



with complex parameters x, y and z. This matrix is invariant under the fi — t symmetry, i.e. the simulta- 
neous exchange of the second and third rows and columns, and under the so-called magic symmetry, for 
which the 22 and 23 entries are written in these particular combinations of x, y and z. 

The second FD mixing pattern that we present is the Bimaximal (BM) scheme 1 18-21 1. Historically 
the BM scenario for lepton mixing was proposed before the TB one, but not in the context of discrete 
non-Abelian flavour symmetry and for this reason in our review comes as second. It is defined by 







V2 V2 

u bm = j 2 2 I ' ( 7 ) 

2 2 V2 

corresponding to the following mixing angles 

sin^ 3 = i, sin0? 2 = i, sin0? 3 = O. (8) 
In this case, the most general mass matrix of the BM type is given by 

mf M = I // r : -; 1 . ,9) 



with complex parameters x, y and z. This mass matrix is still invariant under the n — t symmetry and 
under an equivalent version of the magic symmetry, that defines the 22 and 23 entries. 

Both TB and BM mixing predicts at leading order (LO) a vanishing reactor angle. However, the new 
data on the reactor angle have motivated the analytical research of new mixing pattern that could describe, 
at LO, lepton mixings with a non-vanishing reactor angle. The most complete analysis of this kind have 
been performed in |73 74). However S4 may not been used to predict such patterns J74) and in the 




following we will only discuss suitable modifications of the TB and BM patterns. 

There is a substantial historically difference between the TB and the BM mixing patterns: when TB 
scheme was introduced it was in excellent agreement with experimental data. As a result, a model that 
predicts the exact TB patter as the lepton mixing did not need any corrections, from the phenomenological 
point of view. However in any model so far studied, even in those ultraviolate (UV) completed [75 1, sub- 
leading corrections to lepton mixing are always present. On the contrary, the solar angle of the BM pattern 
has never been in agreement with the experimental data fit: the mixing angles are fitted thanks to a relative 
large correction to the solar angle, of order of the Cabibbo angle A, while only minor corrections affect 
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the other two angles. This makes the BM pattern highly interesting to revisit the idea of quark-lepton 
complementarity 1 76 ■ 80) . 



With the new results on the reactor angle, both TB and BM mixing scheme need corrections to fit the 
data. The origin of such corrections is model dependent but we may catalogue them into two big groups 
once we go in the basis in which all the FD structure arise from the neutrinos: they may arise from the 
diagonalization of either the charged lepton mass matrix or the neutrino mass matrix, once the sub-leading 
corrections are taken into account. Clearly in the most generic case the corrections have both origins. 

Suppose that at LO in the basis in which the charged lepton mass matrix Me is diagonal the neutrino 
mass matrix is diagonalized by the FD mixing matrix Ufd- Thus at NLO we said that we may have 
the following scenarios: 

1) Corrections only from the charged lepton sector: 

M e = M e <°) + 5M e , M v = . (10) 

In this way M e M\ is diagonalized by a unitary matrix V e that may be written as 

1 cf 2 9 cf 3 t 

v, = | -ct*e i c \ 3 9 \ . an 

-cUB -c\%9 1 

where c\ 2 , c\ 3 , c 23 are in the most general case complex parameters with absolute value of order 1. 
According to the TB and BM schemes the final lepton mixing matrix is given by 

U lep = VjU FD , (12) 

where FD stays for TB or BM schemes. The corresponding lepton mixing angles for the TB pattern 
are given by 

(smO^TB =-+fte(c| 3 )0, 

(sin^ 2 ) TB =i-|fte(c? a + cf 3 ) e, (13) 
{smd 13 ) TB =-j= \c c 12 -c e 13 \ 9, 
while for the BM case by 

{sm9^ 3 ) BM =-+ne(4 3 )e, 

{sm9 2 12 ) BM =^ + ^=^e(cf 2 + cf 3 )9, (14) 
{sm9 13 ) BM =-^\cf 3 -cf 2 \9. 

By looking at eqs. ( p"3fl4} for any 9 we may wonder which is the probability to have all the three 
angles in the 3<r ranges allowed by the data assuming the absolute values of the c\ 2 , c 23 and c\ 3 are 
distributed as gaussian centered in one -to be consistent with the assumption they are of order one. 
By plotting such probability as function of 9 we get the Success Rate (SR) of the scheme as shown in 
Fig[T] The 9 value that maximizes the success rate, 9$r, may be seen has the most natural value for 9 
to fit the data. By fixing 9 = 9$r we may visualize the possible correlations among the lepton mixing 



angles by means of scatterplots: for example by looking at eqs. ( 13 H 14 1 we see that 9 23 is expected 
to be only mildly correlated with the other two angles while there could be a correlation between the 
solar and the reactor angles. However if we take the absolute value of the parameters cf • according to 
a normal distribution centered in ±1 the correlation almost disappears as shown in Fig. [2] 
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Q 

Fig. 1 SRs for the TB (left panel) and BM (right panel) schemes at NLO for the NH (cyan) and IH (orange) cases 
according to the parametrization done in eqs | 1 3| 14| assuming that the corrections arise only by the charged leptons. The 
SRs have been obtained by discretizing the interval (0.01, 0.25) for 9 and by taking randomly 10 6 sets (cf 2 , c 13 , C23) 
for each bin. The have been generated assuming a random phase and an absolute value normal distributed around 
1 with a — 0.5. The 9 values that maximizes the SRs are for both the hierarchies 0.103 and 0.107 for the TB and BM 
schemes, respectively. 




Fig. 2 On the upper (lower) line the reactor angle versus the solar (atmospheric) angle for the TB, on the left, and the 
BM, on the right, schemes at NLO according to the parametrization done in eqs |13|14| assuming that the corrections 
arise by the charged lepton sector. The scatter plots have been obtained by fixing 9 = 9sr being 9sr the value that 
maximizes the SR according to Fig[T|and by taking the |cf,-| centered in 1 with normal distribution. Horizontal and 
vertical lines bound the updated 3a ranges for the mixing angles, following the most recent global fit (5). 
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2) Corrections only from the neutrino sector: 

M e = AfJ°> , M v = + 5M V . 

In this case, is diagonalized by Ufd, according to 
a diagonal matrix. In this basis, 8M V can be parametrized by 

ei e 2 
5M V = | ei e 3 
£2 £3 



(15) 

C/J^, • a4 0) ■ U FD , where M" ] is 



(16) 



where the diagonal corrections have been absorbed redefining the diagonal entries of — diag(mi , ni2 , 
In this basis the full neutrino mass matrix 



M v = + SM V , 
is diagonalized by a unitary matrix V v , 



Vu = I -c\l , 



'12 u 
1 

c 23 1 



-13' 

=23 < 
1 



with the |cjy| ~ C(l) and |0| < 1, By writing the ej = e; we easily get 

TO2 — m\ m3 — mi 7713 — mi 

The expression for the mixing angles in the TB case are given by 
1 

2V2 . 
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while for the BM case by 
1 



(sin 0f 3 ) bm 
(sin 0\ 2 )bm 
{sin9 13 ) B M 



2 
1 

= 2 

1 



V2 
^e(c^) 



(21) 



'13 



-231 



The corresponding SRs are shown in Fig. [3] 

As previously noticed, we remind that in general both kind of corrections are present and thus some 
correlations are weakened. However in realistic models more constraints are imposed, as for example 
those relative to neutrino mass splittings and charged lepton masses. In this way the models may maintain 
a good level of predictivity. 

As a conclusion, according to our parametrization and the results reported in Figs. [T] and [3] the BM 
scheme with corrections arising from either the charged lepton sector or the neutrino sector is strongly 
disfavoured, while, on the other side, the TB scheme can be characterized by a success rate up to the 13%, 
with a slightly best result when the corrections arise from the charged lepton sector. These results are in 
agreement with the analysis presented in Ref. JHT) . However, in concrete model realizations it is possible 
that the corrections affecting the lepton mixing matrix are not generic and as a result the success rate could 
reach much higher values: an example we could refer to the model in Ref. |82| , that has been deeply 
analyzed in Ref. ]|8TJ . 
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Fig. 3 SRs for the TB (left panel) and BM (right panel) schemes at NLO for the NH (cyan) and IH (orange) cases 
according to the parametrization done in eqs |20|2f| assuming that the corrections arise only by the neutrino mass 
matrix. The 9 values that maximize the SRs are for both the hierarchies 0.1 and 0.15 for the TB and BM case, 
respectively. The SRs have been obtained as in Fig[T] 




Fig. 4 On the upper (lower) line the reactor angle versus the solar(atmospheric) angle for the TB, on the left, and BM, 
on the right, schemes at NLO according to the parametrization done in eqs. l|20}-(|21}, assuming that the corrections to 
the mixing angles arise by the neutrino mass matrix. The scatter plots have been obtained by fixing 6 = 6sr being 
6s r the value that maximizes the corresponding SR according to Fig. [5] and by taking the |c^ | centered in 1 with 
normal distribution. Horizontal and vertical lines bound the updated 3<r ranges for the mixing angles, following the 
most recent global fit (5j. 
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3 The Group S 4 

The group S4 is the permutation group of four distinct objects, isomorphic to the group O which is the 
symmetry group of a regular octahedron. It may be thought as generated by two generators S, T, satisfying 



S 4 = T 3 = E, ST 2 S = T , (22) 

being E the neutral element of the group. 

It has 4! = 24 distinct elements filled in five conjugate classes according to 

Ci : E; 

C 2 : S 2 ,TS 2 T 2 7 S 2 TS 2 T 2 ; 

C 3 : T,T 2 ,S 2 T,S 2 T 2 ,STST 2 1 STS 1 S 2 TS 2 S 3 TS; (23) 

C 4 : ST 2 : T 2 S,TST,TSTS 2 ,STS 2 ,S 2 TS; 

C 5 : S, TST 2 , ST, TS*, S 3 , S 3 T 2 . 

As consequence there are five irreducible representations that can be read from the character table reported 
in Tab. [2] two one-dimensional and denoted as li and I2, one two-dimensional labelled as 2 and two 
three-dimensional written as 3i and 32- 

Table 2 Character table for 54. \ r Me the five representations and d the five conjugacy classes, n the 
number of elements for each class and h the smallest value for which (% r ) = 1. 



Class 


n 
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x 11 


x l 2 x 2 


x 31 


x 32 


Ci 
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1 
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1 2 
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3 


C-2 


3 


2 


1 


1 2 


-1 


-1 


c 3 


8 


3 


1 


1 -1 










6 


2 


1 


-1 


1 


-1 


c 5 


6 


4 


1 


-1 


-1 


1 



The multiplication rules are given by 

l 1 xR = Rxl 1 =R where R stands for any representation 
I2 X I2 = li 

1 2 x 2 = 2 
1 2 x 3i = 3 2 
1 2 x 3 2 = 3i 

2x2 = li + l 2 + 2 
2x3i=3i+3 2 
2 x 3 2 = 3i + 3 2 



3i x 3i = 3 2 x 3 2 = li + 2 + 3i + 3 2 
3i x 3 2 = l a + 2 + 3i + 3 2 . 

A possible choice for the generators S and T is given by 





(25) 
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where X r — r(X) being r the representation of the generator X and where we omitted the double index 
for T since the representation is identical for both the singlets (triplets). 



With the previous choice of the generators is easy to extract the Clebsch-Gordan (CG) coefficients. In 
the following we use m to indicate the elements of the first representation of the product and (3i to indicate 
those of the second one. The multiplication rules for the 1-dimensional representations are given by: 



11 ® f] = r/ <g> li = 77 with r\ any representation 

1 2 ® I2 = 

l 2 ®3i = 



I2 <E> 3 2 = 




(26) 



The multiplication rules with the 2-dimensional representation are the following: 



2® 2 



2®3i= 3i®3 2 



2 <g> 3 2 = 3i 8 3 2 



with < 



with < 



with < 



11 ~ ai/81 + «2/?2 

12 ~ — ai/82 + a2y8l 

ai/?2 + «2/?i 
ai/3i - a 2 ^ 2 

-|(\/3a 1 /3 2 + a 2 /3 2 ) 
KVWs-c^ft) 

|(\/3Q!2/82 -0:1/82) 
-| (v^a^As + 01^3) 

|(\/3a 2( 82-o 1( 8 2 ) 
-i(\/3a 2 /33 + ai/33) 

-§(\/3ai^2 + a 2 l3 2 ) 
^(V3ai/3 3 - a 2 /3 3 ) 



(27) 
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The multiplication rules with the 3 -dimensional representations are the following: 



3i®3i= 3 2 <g)3 2 = li©2©3i©3 2 



with 



3 2 = I2 



2©3i 



with 



li 



3i 

3 2 
h 



3i 



a x Pi + a 2 f3 2 + a 3 f3 3 

( 7l( a 2^2 - Ck3/33) 

I ^g(-2ai/3i + a 2 /3 2 + a 3 /3 3 ) 



02/33 


+ a 3 /3 2 




+ 0*30! 


«i/3 2 


+ 


a3/32 


- a 2 (3 3 


a 1^3 


- a 3 Pi 


(X2P1 


~ OL X fi 2 



ol\P\ + a 2 (3 2 + a 3 (3 3 

( J g (2a 1 /3 1 -a 2 /? 2 -a 3 /3 3 ) 
I -j^(a 2 (3 2 - a 3 (3 3 ) 




(28) 



We will indicate this basis as the Z?-basis since in this basis both triplets and doublets contract di- 
agonally to give the singlet according to the multiplication rules shown above. Other two S4 basis are 
particularly important in flavour model building and we will refer to them as the TB-basis and the BM- 
basis respectively. TB-basis is characterized by having the T generator diagonal in the triplet and doublet 
representation. All the S4 elements in the TB-basis are obtained by the Z?-basis through a rotation 

X TB = V-X D -V\ 



where Xtb,d ar e me S4 element in the two different basis and V 
for the doublet and triplet representation, respectively. They are given by 



tt/4 



-( 1 

M 1 



—I 

1 




1 1 

■jj UJ 



(29) 

stays for transformations 
(30) 



with uj = e 2 "/ 3 . 

The BM-basis is characterized by having the S generator diagonal in all the representations and is 



obtained by the £>-basis through a rotation equivalent to eq. ( 29 1 

X BM = V-X d -V\ 
where now V = V n , U v jn with 

Utv/4 






i 
-1 



(31) 



(32) 
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In the following sections we will see that there is a tight correspondence between the S4 basis chosen 
and the flavour model building realization. Even if physics is independent by the chosen basis, a framework 
may be easier implemented with a basis instead of with another. 

All the models based on S4 may be distinguished into two major classes: those ones that predict at 
LO a FD mass matrix and thus a lepton mixing matrix that is given by pure numbers, and those ones 
that do not. In the first class we recognize three subgroups: i) models in which the lepton mixing arise 
by the correct combination between neutrino and charged lepton mixing matrix, ii) models in which the 
TB scheme arises only by the neutrino sector, Hi) models in which the BM pattern arises only by the 
neutrino sector. Typically models belonging to subclass i) are realized within the D-basis, while models 
belonging to subclass ii) within the T£>-basis and those of subclass Hi) within the £?M-basis. There is 
an additional class of models that can be thought between these two and it is characterized by models that 
predict the TM mixing pattern. Strictly speaking this pattern is a sort of partial FD mixing because only 
one eigenvector-and correspondingly one mixing angle-is fixed by the flavour symmetry. 

Before describing more in details the models based on the S4 discrete symmetry, we comment on the 
claim [83] that, in order to obtain the TB mixing "without fine-tuning", the finite group must be S4 or a 
larger group containing S4. For us this claim is not well grounded being based on an abstract mathematical 
criterium for a natural model (see also Ref. j84)). For us a model is natural if the interesting results 
are obtained from the most general lagrangian compatible with the stated symmetry and the specified 
representation content for the flavons. For example, in Ref. [85 86], a natural (in our sense) model for the 
TB mixing is built with A4 (which is a subgroup of S4) by simply not including symmetry breaking flavons 
transforming like the 1' and the 1" representations of A 4 . This limitation on the transformation properties 
of the flavons is not allowed by the rules specified in Ref. [83], which demands that the symmetry breaking 
is induced by all possible kinds of flavons (note that, according to this criterium, the SM of electroweak 
interactions would not be natural because only Higgs doublets are introduced!). Rather, for naturalness we 
also require that additional physical properties like the VEV alignment or the hierarchy of charged lepton 
masses also follow from the assumed symmetry and are not obtained by fine-tuning parameters: for this 
actually A4 can be more effective than S4 because it possesses three different singlet representations 1, 1' 
and 1". 



4 5 4 -Based Neutrino Models 

In this core section we give an up to date overview of the S4 models present in the literature according to 
the classification we did at the end of Sec. [3] between models that predict an exact FD lepton mixing at 
LO, and those that do not. We already said that in the first group there are the TB and BM mixing models, 
while the TM mixing one is in the middle and we include it into the second group for reason that will be 
clear in the next sections. 



4.1 FD Models 



Here we introduce the theoretical motivation behind the FD mixing patterns as already discussed in 1 73 



74 83 87 88 1. The same arguments apply not only to S4 but to any non-Abelian discrete symmetry. For 
this reason we will be as generic as possible following 1 73 , 74 1 . 

Let us consider a non-Abelian discrete group G/, of elements {g^} with i = 1, da , being do the 
dimension of the group. Let r be an irreducible representation of dimension n. The n x n matrices r(gt) 
will be the representations of the elements of Gf. Suppose that Gj has a given set of subgroups Gf t , 
I = 1, .., h of dimensions dc r Then each subgroup is generated by a set of do, {gi k } subset of {gi}. A 
n x n matrix M is invariant under the subgroup G/, if for each gi k , k = 1, .., da r we have that 

r(g lk ) T ■ M ■ r(g lk ) = M , (33) 

if M is symmetric, or 

r( gi Jt . M ■ r(g, k ) = M , (34) 
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if M is hermitian. In general the irreducible representation r of Gf is a reducible representation of the 
subgroup Gf, and decomposes into ri 1 © • • • © ri m irreducible representations of Gt v This means that 
exists a unitary transformation SI; that rotates each r(gi k ) into a block diagonalizable form corresponding 
to ri 1 (gi k )©•••© n m (gi k ), or in other words that 

a\-r(g lk )-n l =f(g lk ), (35) 
where we have indicated f(gi k ) = (gi k )©•••© r\ m (gi k ). Thus eq. ( 36 1 becomes 

f(g lk ) T -M-f(gi k ) = M, (36) 
where we have indicated 

M = flf ■ M ■ Qi . (37) 
Considering eq. (|34|> we straightforwardly get 



f(g lk f ■ M ■ f{g lk ) — M , M = Sl\ ■ M ■ . (38) 

M maybe interpreted as the matrix M in the G /, basis, being G /, the subgroup of Gf under which it is 
invariant. 

Let us consider now the basis in which both neutrino and charged lepton mass matrices are diagonal. 
Since no degeneracy is present in any of the two sectors, if there is any residual flavour symmetry this has 
to be Abelian. Assuming that light neutrinos are Majorana particles such Abelian symmetry is limited to 
be Z2 x Z2. The reason is very simple: the Majorana nature of neutrinos allows them to be charged only 
under a Z 2 , giving them and odd or even flavour parity. However with only one Z 2 , it is not possible to 
have a diagonal neutrino mass matrix -we are in the mass eigenstate basis- and therefore it is necessary to 
introduce a further Z2 to forbid off -diagonal entries. For what concerns the charged leptons, being Dirac 
particles, there is a infinite class of Abelian symmetries, both continuos and discrete, that may be chosen 
to get a diagonal charged lepton mass matrix. The minimal choices are obviously Z2 X Z2, Z2 x Z2 x Z2 
and Z3. 

Coming back to the general discussion, it is clear that G / has to be read as the main flavour group, while 
the subgroup G f t are only its Abelian subgroups that define the flavour structure of neutrinos and charged 
leptons. As a consequence the representation r; . (gi k ) are all one dimensional and consequently both r{gi k ) 
and M are diagonal. Finally S7; corresponds to the matrix that diagonalizes M. In other words if M is the 
neutrino (charged lepton) mass matrix M v (MJM e ) invariant under Gf v (G/ e ) subgroup of Gf, the lepton 
mixing is given by 

U lep = n\ ■ SI, . (39) 

We automatically get the important information that G e cannot coincide with G v because otherwise Ui ep 
would coincide with the identity matrix. 

Any FD lepton mixing is obtained by choosing a non- Abelian discrete G / and consequently its sub- 
groups G„ and G e for the neutrino and charged lepton sector respectively, that fixes unambiguously SI, 
and Sl e . 

Finally a remark: when we say that we fix a subgroup Gf t of Gf generated by the subset {gi k } we 
mean up to subgroups generated by elements belonging to the same classes of {gi k }. Since all the elements 
belonging to the same class may be obtained by acting on only one element of that class, it turns out that 
the subgroups originated by elements of the same class give rise to a series of SI; that are identical up to 
permutation of raws and columns. 

Let us now consider the case of S4. by looking at Tab. [2] and according to what said so far, the residual 
symmetry G v ~ Z2 x Z2 in the neutrino sector has to arise by one element of C2 and one element of 
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C4. If we take S 2 for the first and TST for the latter and impose that M v is invariant under both of them, 



according to eq. ( 36 1 we get 




M v = b c . (40) 



This matrix is diagonalized by a maximal rotation in the sector 23 by means of R23 (tt/4) using the notation 
introduced in eq. pj. According to what we have said, taking other two elements of the same classes we 



expect the same structure of eq. (40 » but involving different columns/raws: indeed if we consider TS 2 T 2 
of C 2 and STS 2 of C 4 , M v is diagonalized by R 13 (n/4), while, taking S 2 TS 2 T 2 and ST 2 , the rotation is 
Ri 2 (-k/A). Finally all the other combination lead to the same structures or to the limiting cases for which 



the off -diagonal entries are put to zero, that is c = in eq. (40 1. 

For what concerns the charged lepton sector we have now two options: we may take the subgroup 
generated by an element of the C3, corresponding to G e ~ Z3, or by an element of C5 corresponding to 
G e ~ Z 4 . By taking T of C 3 and requiring that M e M\ is invariant under T, where M e is the charged 
lepton mass matrix in the LR notation, forces 



M e Ml = \ g* f g , (41) 




that translates into 

(h hi h 2 
h 2 h h x I • V R , (42) 
hi h 2 h 

being Vr an arbitrary non-physical rotation in the right handed charged secto^] The latter is diagonalized 
by the matrix ~U W introduced in Sec. [3] to switch basis from the D-basis to the TB-basis and it has been 



defined in eq. (30 1. It is clearly not surprising that in this case by combining U u with R 2 3(tt/A) according 
to 

U lep = Ul-R 23 (n/4), (43) 

the final lepton mixing is the TB one. Notice that in the TB-basis the T generator element of C3 is diagonal, 
thus M e is diagonal and all the mixing pattern arises from the neutrino sector. From this the name of the 
basis. 

Let us now consider the case in which G„ is generated by the elements S 2 TS 2 T 2 of C 2 and ST 2 of C4. 
We already said that this choice leads the following neutrino mass matrix 

a b 

M v = I b a I , (44) 

0c 

diagonalized by Sl„ ~ i? 12 (7r/4). In this case, for the charged leptons, we can take G e ~ Z 4 : this means 
that we require M e AfJ to be invariant under, for example, the action of S, element of C5. This leads to 



M e Mt = g ih , (45) 




^ Notice that Vr is non-physical as long as we deal with the SM. In a supersymmetric context Vr may be related to the corre- 
sponding rotation in the right handed charged slepton sector with no trivial consequences for the phenomenology of the model. 
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where /, g and h are real parameters. This matrix is diagonalized by a i?23(— 7r /4) up to a phases: more 
precisely by: 



1 











1 


1 




f 2 





V2 


V2 










)■: 


— i 


• 










i? 23 (-7r/4) = u 75 -7i • -z . (46) 



In this way by combining eq. (46 1 with i?i 2 (7r/4) arising by the neutrino sector we get the BM pattern of 



eq. |7]). The diagonal phase matrix may be absorbed by a field redefinition. 



Let us concentrate on the consequence of eq. (45 1 for M e once we have absorbed the phases of eq. ( |46) . 
M e presents the generic structure 

/ goe 1 ^ \ 

M e = hoe** h 2 e^ 2 • V R , (47) 

\ hie i( - mir+ ^ y/h\ + hl- hje^+r*) J 

where g , h , hi , h% are real parameters, ip g i their phases, m,n integer numbers and Vr an arbitrary 
rotation in the right handed sector. Here we have considered the case h\ + — h\ > 0. However the case 
h\ -\-h\ — h\ < is easily got by changing the sign and replacing n with n + 1. M e describes the charged 
lepton mass matrix, so we need 3 different eigenvalues. For the most generic case M e MJ has three distint 
eigenvalues 

2 2 

m 1 =g , 



m% 



2 = hi + hi - y h\ + hlh\ + h\{Kl - h\) + 2h 2 h 1 h y h\ + h% - h\ cos(to - n)n , ^g-j 

ml = hi + hl + J}4 + h\h\ + h\{h\ - h\) + 2h 2 h 1 h ^ h\ + h% - h\ cos(m - n)n . 

However in many models, it is easier getting high symmetric mass matrices: consider for example the case 
in which ho = hi. In this case the two last eigenvalues becomes 



777.2.3 = h 2 + h - (±) V h 2 + h + 2h 2 h COs(?77 - 77)tT . (49) 

In order to have ml ^ Owe need m — n = 1 + 2kir, with k an integer. Thus taking m = 1, n = and 
k = -and Vr = 1- leads us to the simplest form for M e we may get in the D-basis 



M e = h h 2 , (50) 




where now go , hi are complex. This texture for the charged lepton mass matrix has been used in Ref. 
Notice that in that context it was deduced by simple phenomenological motivations. On the contrary in the 
D-basis the most generic charged lepton mass matrix invariant under Z 4 is not easily obtained in a concrete 
model. Once again the choice of the basis helps: in the BM-basis S, generator of Z4 is diagonal, thus M e 
is diagonal with 3 distinct mass eigenvalues and all the lepton mixing arises by the neutrino sector. 

Recently, on the basis of only theoretical group theory considerations, the authors of Ref. ]90| showed 
that in the case of 54 new patterns predicting 6*13 7^ may be obtained. This could seem in contrast with 
what we have showed so far and with what is stated in Refs. |73 74|. However, this is not the case. Indeed, 
under the assumption that neutrinos are Majorana particles, their mass matrix presents always a Z 2 x Z 2 
symmetry. One of the two Z 2 could be accidental and not belonging to the flavour non-Abelian group. It 
has been shown in Ref. |90|, from general group theoretical considerations, that certain correlations among 
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the mixing angles follow from the residual symmetries of the charged lepton and neutrino sectors, once the 
initial flavour symmetry is broken. As a result, the relations showed so far and in Refs. [73 74) are special 
cases of the more general expressions recovered in Ref. [90]. Although we find the results in Ref. [ 90 1 
interesting and potentially inspiring for new strategies in neutrino model building, the construction of a 
concrete model predicting the correlations showed in Ref. [90| could be quite challenging: indeed, when 
only one of the two Z 2 symmetry factors of the neutrino mass matrix is contained in the flavour group 
Gf, it is necessary the accidental presence of a Z 2 factor aside the flavour group, Gf X Z%, and this is 
subjected to a well precise alignment in the flavour space of the flavon vevs and to a specific choice of the 
transformation properties of all the fields under Gf. 



4.2 Breaking the Flavour Symmetry 

From a group theory point of view, it is quite natural to assume that a given non-Abelian discrete group G / 
breaks down to two different subgroups G v and G e in the neutrino and charged lepton sector, respectively. 
However from a practical point of view it is not easy that a model accounts for the correct mechanism to 
break the symmetry. 

Typically the breaking of the flavour symmetry is felt by the SM fermions, with the addition of the right 
handed neutrinos, through their couplings with scalar fields that transform not trivially under the flavour 
symmetry, usually addressed as flavons. In the most traditional scenarios the flavons are singlets of the 
SM gauge groups and bring only flavour charges, while the SM boson Higgs field is not charged under 
the flavour group, or better said it is not charged under the non-Abelian part of the flavour group, since 
it is often charged under additional Abelian symmetries. Additional scalars necessary to induce neutrino 
masses (like SU(2) l triplets) behave in general as the Higgs field with respect to the flavour group. 

Consequently, these models are written in terms of non-renormalisable effective operators suppressed 
by suitable powers of the cut-off. Since in general the higher-dimensional operators cannot be neglected, 
it is necessary to take into account at least the NLO corrections -sometimes even NNLO- as long as an 
UV renormalizable version of the models is not studied. In the most traditional spirit we expect that the 
flavons develop VEVs inducing the spontaneous breaking of the original flavour symmetry, thus giving rise 
to fermion mass textures. However one can think to ED scenarios in which the flavons leave in the bulk 
while the SM fermions on the boundaries. In this case adequate boundary conditions (BCs) could break 
the symmetry and make the desired job. In this brief discussion we will not consider this latter case, but 
we acknowledge it as an alternative to the traditional dynamics of the flavour symmetry breaking, that in 
our opinion is one of the weakest point of these models. 

Let us consider the spontaneous breaking of the flavour symmetry: since G / needs to be broken down 
to G v and G e in the neutrino and charged lepton sector, respectively, and G v ^ G e , we need two distinct 
set of flavons, denoted as {<fi u } and {0 e }, that implement the breaking chains Gf G v and Gt — > G e 
and communicate the breaking of the symmetry to neutrinos and charged leptons. We used the notation 
{4>v.e} to be as general as possible, meaning that for each breaking we may have more than one flavon 
field. If now, for concreteness, we consider a single flavon <f>i ~ r, where r is an irreducible representation 
of Gf, denoting as (<f>i) its VEV, we have that Gf —> Gf t if for any r(gi k ) it holds 

(4>i) • r(gt k ) = r(gi k ) , (51) 

being {gi k } the elements of Gt r For example, if we work in the Z?-basis a triplet (j> of S4 induces the 
breaking S4 — s- Z 2 x Z 2 where Z 2 x Z 2 is generated by S 2 and TST if {4n) is aligned as (1,0,0). 
Similarly cj) induces the breaking S4 —> Z3 generated by T if its VEV aligns as (1,1,1). 

Given a set of flavons, the spontaneous breaking of Gf to one of its subgroup is naturally realized: 
typically we need a negative mass square term in the flavon scalar potential and then the relative size of the 
potential parameters determine the subgroup in which G r breaks. On the contrary, in the most general case, 
Gf does not break to two distinct subgroups: this turns into the necessity to separate the scalar potential 
of {4>u} from that of {<fi e }- The easiest way to realize this is by introducing additional Abelian symmetries 
that also prevent {</>„} coupling to the charged lepton sector and {4> e } coupling to the neutrino one. A 
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further help is provided by promoting the models to be supersymmetric, even if it is not mandatory [91 ]: 
in the supersymmetric case by imposing the U(1)r continuous symmetry and by introducing an additional 
set of scalar fields, usually called driving fields, the desired breaking chains are relatively simply realized. 

The drawback of splitting up the potential is that it generates accidental continuous symmetries that 
once broken give rise to massless goldstone bosons. The problem of the massless particles is avoided by 
introducing soft flavour symmetry breaking terms, but the origin of such soft terms has not been addressed 
so far. 

Finally we remind that the scalar potential may be splitted embedding the models in large EDs: in this 
way the two sets of flavons live one on the IR brane and the other on the UV one and the interference 
between the two sectors are suppressed by the size of the extra dimensions. Notice that this approach is 
completely different with respect to breaking the flavour symmetry by means of BCs: indeed in this case 
the EDs setting is used to separate the two set {cf) u } and {<p e }, making them living into different branes 
and the flavour symmetry breaking is induced by flavon VEVs as usual. 

In conclusion the correct flavour symmetry breaking necessary to induce FD mixing patterns require 
additional ingredients that make the models more baroque and less appealing: extra symmetries, extra 
dimensions, extra scalar fields, unexplained soft terms. This is a serious drawback for FD mixing models, 
because as much as natural they seem under group theoretical considerations as less natural they look after 
a concrete model building realization. For these reasons no FD mixing models represent a complementary 
and compelling approach to account for the flavour puzzle: the underlying flavour symmetry still maintain 
a certain level of predictivity but its breaking is easily realized. 



4.3 General Criteria to Classify Flavour Models 

In the next sections we will summarize the main features of the models that at LO predict TB and BM 
lepton mixing, respectively. In order to discriminate among them and to identify the most complete ones, 
we list a series of general criteria typically addressed in flavour model building: 



a) realistic realization: we have said that flavour breaking patterns are not easily realized. In principle 
one could build a flavour model assuming unjustified vacuum alignments hoping that it will exist a 
scalar potential or a mechanism to explain them. Models that do not enter into the details of flavour 
symmetry breaking are incomplete and weak because based on a hypothesis that could not be true; 



b) predictions in the neutrino sector: there is a plethora of way to explain the origin of neutrino masses 
that reflects in different phenomenologies. Moreover the flavour symmetry may be implemented in 
such a way to get interesting predictions for neutrino observables, such as 0^2/3 decay, the kind of 
hierarchy, the reactor angle and CP-phases. Predictions are a plus value for this kind of models; 



c) GUT embedding: it is quite appealing to try to combine flavour symmetry with grand unified theories 
(GUT). This is deeply challenging because in GUT models quarks and leptons are (partially) unified 
while quark and lepton mixing and mass hierarchies are completely different. Models that give a 
unified description of quarks and leptons can be considered more complete. 



d) quark sector: being less ambitious than at point c), we may wonder if quark mixing and mass hierar- 
chies may be described by the same flavour symmetry as leptons, even if not in a GUT context; 



e) further predictions: models satisfying criteria c) or the simpler d) may have interesting prediction in 
the quark sector (such as the Gatto-Sartori-Tonin (GST) | |92) or the Georgi-Jarlskog (GJ) (93) rela- 
tions) or interesting correlations between charged lepton and down quark masses. Moreover fermion 
masses may satisfy sum rules. 
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Table 3 Model properties according to the criteria enunciated in Sec. |4,3| FD models are grouped into two 
classes, namely TB and BM, referring to those predicting TB and BM patterns as the LO lepton mixing 
matrix, respectively. In a third group, no FD, we put all the S4 flavour models that do not have at LO a FD 
lepton mixing matrix. " (0, h)" stays for the flavour and Higgs scalar potential study; "u" indicates the 
neutrino sector properties; "GUT" if the model has a GUT framework; "5" if it describes also quark sector; 
"Qrei" if it predicts quarks relations such as GST and/or GJ; finally under "extra" we grouped all the other 
relevant features. If the latters are absent we put a For what concerns neutrino porperties "type-X" stays 
for the kind of See-Saw, "ISS" for Inverse See-Saw and "W.op." for Weinberg operator. 





v{4>, h) 


V 


GUT 


q q re i 


extra 


TB 



|94 

m 

|55_ 
5 _+_ 
98 



97} 



rorji 

[102| 
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roughly 
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no 
yes 
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no 
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BCs in ED 
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type I 
type I+II 

type I 
type I+II+III+W.op. 

type I 
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no 
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SU(3) embedding 


no 
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no 


complicated scalar sector 
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no 




no 
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yes 
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8D,orbifolding 


5/7(5) 
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5*7(5) 
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no 


only Cabibbo in Vqkm 
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yes 


no 




no 


no 


no 
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1 88 yes lower bound for mi no no no 

55 yes type II PS yes yes QLC 

[ lOo I no type II 50(10) yes no QLC, renormalizable 

[107 1 yes W.op. SU(5) yes no b — r, 5D orbifolding 



No FD 



11081 
TO 

(no I 



yes 


type I 


no 


yes 


no 


embed, in 50(10) x G f 


yes 


W.op., e atm = tt/4 


no 


yes 


no 


tx product of many Z 2 


yes 


type I, TM 


no 


no 


no 


new sum rule 


roughly 


ISS 


no 


yes 


no 


new sum rule 



4.3.1 TB Models 



Here we deal with the models that predict TB mixing at LO order. Armed by the criteria exposed in Sec. 4.3 
and by Tab|3]it is easier to describe the different features and properties. 

In all the models the flavour symmetry is represented by a product of 54 and additional Abelian sym- 
metries: only in the model in Ref. |103| , a Frogatt-Nielsen continuous symmetry is introduced, while in 
all the others discrete symmetries have been used. Moreover, all of them, a part from the realisation in 
Ref. fl9T) , are supersymmetric. As a result, we did not list these features in Tab|3] 
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The models in Refs. |95 1 and [99] do not satisfy criterium a), while in Ref. [94 1 only a rough analysis has 



been done. In this paper an embedding in SU(3) is necessary to get the so-called Koide relation 1 1 12 - 114) 
among the charged lepton masses, 

m e + rrifj, + m T _ 2 

- - 3 , P^J 



m e + v m ,, + v m T 

that is very well verified for pole masses. However, the VEV alignment is ad hoc. In the model in Ref. [95 1 
neither the flavour nor the GUT scalar potentials are studied and the scalar sector to get fermion charged 
masses is deeply cumbersome. Similar situation for the model in Ref. [99]: the scalar potentials are not 
studied and neutrino masses are assumed to be generated by the type-II See-Saw mechanism. On the 
contrary, for all the other models the scalar potential is analyzed in details. In particular, in the model in 
Ref. 1 102 1, a really exhaustive study of both the scalar potentials is performed, both at LO and NLO (in 
Ref. H 1 1 1 1 a variation of the model in Ref. |102[ has been presented with the focus on the CP violating 



Dirac phase). In Ref. |101| , S4 is obtained by orbifolding a 8-dimensional space and then the flavour 
symmetry is broken by BCs. 

For what concerns point b), in the majority of the models, neutrino masses are induced by type I See-Saw 
mechanism. However, in Ref. [99], the type II mechanism is implemented, while in Ref. p04[ , there is an 
interplay between type-I and type-Ill. In Ref. J9T[ , neutrino masses are induced by the interplay between 
type-I and II See-Saw. Refs. J%) - (97) collect different scenario those more mechanisms to generate 
neutrino masses are involved. All the models have NLO corrections that may give a #13 in agreement with 
the most recent data. However, NLO corrections are under control only in those models that satisfy criteria 
a). Quite interestingly in Ref. [ 105 ], only the NH case with mi = is predicted. 

All the models in Refs. ]95|99|10T[|104 ) are embedded in a GUT scenario: the first four in SU (5), while 
the last two in S'O(IO). All the models fit quark masses, even if the model in Ref. | |95) has a complicated set 
of scalars to fit the charged masses. For what concerns the mixings, the model in Ref. 1 103] predicts only 
the Cabibbo angle, while the other models fit all the CKM matrix. However, only the model in Ref. 1 101) 
satisfies the GST and GJ relations, thanks to a specific set of messengers. 



4.3.2 BM Models 



The number of models that at LO predicts exact BM mixing in the context of S4 is relatively small. All the 
models are embedded in a supersymmetric scenario and have additional Abelian symmetries: we did not 
put these features in Tab. [3] as we did not for the TB mixing models. As done in Sec. 4.3.1 



we summarize 



them by means of the criteria given in Sec. 4.3 



In Ref. 1 106 1, a renormalizable model with the addition of messenger fields is built. However, both 
the flavon and the Higgs scalar potential are not studied. On the other hand, all the other models in this 
category satisfy criterium a). 

For what concerns the neutrino sector, in the model in Ref. 1 88 1 a lower bound for the lightest neutrino 
is recovered, in the context of the type I See-Saw. In the models in Refs. [89) and fj06|, the dominant 
contribution to neutrino masses arises from the type-II See-Saw: both the models are embedded into a 
GUT context, Pati-Salam and 50(10), respectively, and therefore in general both type-I and type-II See- 
Saw contributions are present; the dominance of type-II must be suitably justified. However, only the model 
in Ref. [89 1 studies in details the Higgs scalar potential to check the validity of the type-II dominance, while 
it is simply assumed in Ref. 1 106 1. Finally, in Ref. [107], neutrino masses can be induced by either the 
Weinberg operator or the type-I See-Saw. 

In addition to the models in Refs. [89| and 1 1 06 1 , also that one in Ref. |107| is embedded in a GUT 
scenario, that in this case is SU(5). Moreover this model presents an extra dimensional setup with the fifth 
dimension compactified on an orbifold, mainly to provide an explanation to the doublet-triplet splitting 
problem. 
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The non-GUT model, in Ref. [88 1, does not deal with quarks, while for the other models, the quark 
sector is automatically accounted for by the GUT embedding. In particular, the model in Ref. [89] presents 
at LO the GJ relation at the GUT scale, slightly broken once NLO effects are introduced. In Ref. [106 1, 
the GJ relation is reproduced only approximately even at LO, since m M = — 3m s is obtained assuming the 
smallness of a contribution. In Ref. |107[ , only the b — r unification is obtained at the GUT scale but it is 
preserved even at NLO. Finally in all the models, a part for that one in Ref. (106) , a Froggatt-Nielsen U{\) 
symmetry is introduced to explain the charged fermion mass hierarchies. 



4.3.3 No FD Models 



As already stated in Sec. 4.2 flavour models based on a discrete non-Abelian flavour symmetry may 



maintain a good level of predictivity even if they do not reproduce a FD lepton mixing at LO. This is due 



to the reduced number of free parameters dictated by the flavour invariants. All the criteria in Sec. 4.3 may 
be extended even at these models, even if point a) is clearly more easily realized. 



In Ref. (108), it is presented a model embeddable in a GUT 5*0(10) x G f , being G f ~ SO(3) f or 
SU(3) f continuous flavour symmetry. The low energy Higgs doublets transform not trivially under £4, 
thus the flavon and the Higgs sector need to be treated together and the S4 breaking scale coincides with 
the EW one. Two numerical examples are provided for which the mass matrix textures lead to realistic 
lepton and quark masses and mixings. Furthermore, it is discussed the Higgs potential that allows for the 
appropriate VEV configurations required for such textures. 

The model in Ref. [ 109] is a supersymmetric model in which the flavon and the Higgs sectors are treated 
together: a large set of SU(2) Higgs doublets transforming non-trivially under the flavour symmetries is 
present. Indeed the model have six Higgs doublets and the minimization of the potential gives a residual 
fi — t symmetry that predicts maximal atmospheric angle and vanishing reactor angle. Both quarks and 
leptons mass hierarchies and mixings are fitted, and neutrino masses are induced by mean of the Weinberg 
operator. However, the full flavour symmetry of the model is baroque, being Gf = £4 x Z$ q x Z2 q X 
{Zi2e x x Zir) and represents a weak point of the model. 

In Ref. p8[ , it is described an S4 supersymmetric model where the scalar VEVs break S4 —> Z2 in 
the neutrino sector and to Z3 in the charged lepton one. In this way, when charged lepton are diagonal, 
an exact TM lepton mixing is obtained and the mixing angles are expressed in term of deviations from 
the TB ones, namely s, a and r for the solar, atmospheric and reactor angle respectively. The interesting 
smoking-gun signature of such scenario is s ~ and the sum rule 2a + r cos S ~ 0, being 5 the CP-Dirac 
phase. 

Finally Ref. |1 10[ presents a supersymmetric model that justifies from first principles an interesting 
new sum rule for neutrino masses in the context of the inverse See-Saw (ISS) mechanism. Once more the 
drawback of the model is a very large flavon sector and consequently a really complicated scalar flavon 
potential. 



5 Phenomenology 

In this section we deal with the phenomenological signatures of the models described in the previous 
sections. Without entering into details of a specific model, we present a general analysis that could be 
used as a guide to perform a more detailed study. We mainly concentrate on the presence of specific sum 
rules for the neutrino masses, on some neutrino observables, such as the sum of the neutrino masses, the 
kinematic electron neutrino mass in the single beta decay and the neutrino-less double beta decay effective 
mass, and on some lepton flavour violating transitions, like /i — > ej. These observables could be very 
useful to test the validity of flavour models and to distinguish one model from the other. 
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5.1 Sum Rules and Related Observables 



In the context of FD mixing pattern, the neutrino mass matrices depend only on three complex parameters 
(cfr. with eqs. (|6]i and (|9]l): three real parameters describe the neutrino masses and the remaining two 
corresponding to the Majorana phases. Indeed, the three mixing angles are fixed and independent on the 
parameters of the mass matrices and the Dirac CP phase is undetermined, being #13 = 0. 

Often, in flavour models based on non-Abelian discrete symmetries, only two independent complex 
parameters enter the neutrino mass matrix and therefore the three neutrino masses undergo to a well- 
defined correlation. Up to now, four different sum rules have been presented in literature 1 1 15 1 16) : 



mi = a m 2 + ft 

1 1 o 1 

= a VP 

m\ 777 2 7B3 

/mi = ol \Jvri2 + P \frnj, 

1 1 a 1 

1= = a—=+ (3 —= 
/mi v m 2 v m 3 



(53) 



where mi are the neutrino mass eigenvalues and a and f3 two real and positive parameters, depending on 
the details of a specific model. 

The first of these relations is usually present in model where the neutrino masses are described by the 
Weinberg Operator or in the context of the type II See-Saw mechanism (see for example Refs. [96,97]); 
the second and third sum rules are instead described in model with type I and type III See-Saw mechanism 
(see for example Refs. [97 ,98 1); while the last relation arises in the context of the inverse See-Saw (see for 
example Ref. (TlOj). 

Once a specific mass ordering is chosen, it is possible to express the heaviest states in terms of the 
lightest ones and the mass squared differences: in the case of the normal ordering (NO), one can write 



m 2 



''sol ' 



m 3 



(54) 



while in the case of inverse ordering (IO), they are 



mi 



777 2 



A m? , 

SOL 



(55) 



As a result, entering these expressions into the sum rule of a specific model, all the three neutrino mass 
eigenvalues are determined. 

The precise determination of the three neutrino masses at experiments would represent a strong con- 
straint for neutrino flavour models, testing for each model the corresponding sum rule. However, the 
present (and the expected future) sensibility is too low to allow such test. On the other hand, the pres- 
ence of these correlations among the neutrino masses have an interesting impact on other observables: the 
sum of the absolute neutrino masses E, the kinematic electron neutrino mass in beta decay 777,3 and the 
neutrino-less double beta decay effective mass (m ee ): 



^2 777 fe , 



k=l 



\ k=l 



E^e 2 



'k m k 



c 12 c 13 7771 + s 12 c 13 e 2 1 777 2 + Si 3 e 31 777 3 I , 



(56) 
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where U e k are the elements of the first row of the lepton mixing matrix, c,j and refer to cosines and 
sines of Oij, while 021.31 are the Majorana phases in the usual convention. 

Taking the two Majorana phases in the range [0, 7r], one can generate the plot in Fig. [5ja), where the 
0^2/3 effective mass is drown as a function of the lightest neutrino mass for both the mass orderings. 
The horizontal lines corresponds to the past and future sensitivities of the Qv2j3 decay experiments: the 
Hidelberg-Moscow experiment |1 17) put an upper bound at 90% of C.L. of (m ee ) < 0.21 — 0.53 eV, 
while a part of the collaboration claimed the observation of the 0^2/3 decay [ 1 1 8 1 corresponding to a value 
(m ee ) = 0.32 ± 0.03 eV; the expected future sensitivities of GERDA II fl!9j(GERDA III, EXO fl20) ) 
are of (m ee ) - 90 - 150 (10) eV, while that one of CUORE (TIT) is of (m ee ) - 41 - 96 (10) eV. The 
vertical ones to the future sensitivity of KATRIN fl22| , 0.2 eV, and MARE fll3) , 0.1 eV. 

In Fig. |5|b), correlations among the sum of the neutrino masses, the single j3 parameter and the 0v2/3 
effective mass are shown. The vertical lines show the recent bounds on E |124| considering two contexts 
and different combination of the cosmological data: the Blue corresponds to the standard cosmological 
model ACDM with massive neutrinos and the Red the generalization with non-vanishing curvature and 
with u> 7^ —1 in the Dark Matter equation of state. CMB stands for the Cosmic Microwave Background, 
HO for the Hubble Constant, SN for the high-redshift Type-I SuperNovae, BAO for the Baryon Acoustic 
Oscillation, LSSPS for Large Scale Structure matter Power Spectrum. 



- GERDA I 
CUORE 



Heidelberg-Moscow 



Klapdor Claim @3(r 



GERDA 111 - EXO 




m, (eV) 
(a) 0v2f3 effective mass. 



S(eV) 
(b) Correlations. 



Fig. 5 (a) The 0^2/3 effective mass as a function of the lightest neutrino mass for both the mass orderings: in Red 
(Blue) the NO (IO). The coloured areas (dashed lines) correspond to the neutrino oscillation parameters at la (3a) 
from Ref. [5 1. (b) Correlation between £ and mp and (m ee ). See details in the text. 



Considering a specific model and the corresponding sum rule, it is possible to restrict the parameter 
space of the three observables: indeed the two Majorana phases enter the sum rule and prevent the deter- 
mination of a single point in their parameter space. Studies on such bounds in specific models have been 
presented in Refs. ]1 10||1 15p 16], where it has been shown how the implementation of the sum rules could 
represent a strong improvement to test flavour models for neutrino masses and mixings. 

However, as also pointed out in Ref. |115| , these correlations among the neutrino masses are perturbed 
due to the presence of the higher order operator contributions. As a result, it is not possible to identify 
a precise point in the parameter space (£, mp) or a precise lower bound on (m ee ), the uncertainty being 
proportional to the amount of the corrections. For the FD mixing patterns discussed in the previous sections 
the required correction to the reactor angle to suitably reproduce the corresponding experimental value 



Copyright line will be provided by the publisher 



fdp header will be provided by the publisher 



23 



must be of the order of O(0. 1). As a result, the relative corrections to the neutrino masses are expected to be 
of the same order of magnitude, that introduce a sensible uncertainty to the sum rule and the corresponding 
constraints on E, mp and (m ee ). 

5.2 Lepton Flavour Violation 

In this section we discuss the implications for lepton flavour violating (LFV) processes of the classes of 
models described in the previous sections. In particular, we will focus on radiative lepton decays, that turn 
out to give very strong constraints. 

For all the effective models described here, we will consider a working framework in which SUSY 
is broken at a high-scale, higher than or comparable to the flavour breaking one Af, in a hidden sector. 
Subsequently, this breaking is transmitted at the lower scale A of few TeV by a mechanism such as gravity 
mediation. In this context, the soft SUSY breaking terms undergo the flavour symmetry and specific non- 
universal BCs arise at the scale Af. It is the non-universality of the BCs that differentiates the contributions 
to the radiative lepton transitions of the flavour models under discussion from the well-known MSUGRA 
framework. 

5.2. 1 The SUSY Lagrangian and the Soft SUSY Breaking Terms 
The Lagrangian of a generic SUSY model is given by 



JSf = J d 2 0d 2 Q K(z, e 2V z) + J d 2 Sw(z) + h.c 



d 2 Of{z) WW + h.c. 



(57) 



where /C(z, z) is the Kahler potential, w(z) is the superpotential and f(z) is the gauge kinetic function. V 
is the Lie-algebra valued vector supermultiplet, describing the gauge fields and their superpartners while 
W is the chiral superfield describing, together with the function f(z), the kinetic terms of gauge bosons 
and their superpartners. 

When the flavour symmetry is broken, the kinetic terms could receive non-canonical contributions. In 
totally general notation, we can write 



&Hn = i K iS li o» tj + i if P. t\o» Dp I) + Kn D» it If + K% DH\ I) , (58) 

where is the covariant derivative and K^> are hermitian 3x3 matrices in the flavour space. 
The superpotential can be written distinguishing several parts: 

w = we + w v + Wd + Wh , (59) 

where, once the flavour symmetry is broken, we (w v ) is responsible for the charged lepton (neutrino) 
masses, Wd contains the self-interactions of the flavons and possibly their interactions with other super- 
fields, apart from those already described in we and w u and with the Higges. Finally, the last term Wh is 
associated to the /x term 

w h = Xfj, H u H d . (60) 

The soft SUSY breaking terms are generated from this SUSY Lagrangian analytically continuing all the 
coupling constants (such as the couplings in the superpotential and in the Kahler potential) into superspace, 



i.e. by promoting all the coupling constants to superfields with constant 2 and 9 2 8 2 components 1 125 
126) . In particular, the Lagrangian describing the slepton masses is 



J2f m D (e e c )M 2 e [ t )+vm 2 vLL v, (61) 
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where 



m eLL m eLR 



Mi = t ulJ e ^ n , (62) 

V m eRL m eRR J 

with m 2 e V \ LL and m 2 eRR being 3x3 hermitian matrices, while m\ LR = (Tn 2 RL ) . Each of these blocks 
receives contributions from different part of the SUSY Lagrangian: 

™Zrr = {™Irr)k + (m 2 eRR ) F + (ml RR ) D , (63) 

m lRL = ( m lRL)l + (™eAE,)2 ■ 

The indexes specify the origin of these terms: K stands for the Kabler potential; F and D for the SUSY 
F— and D— term contributions; (m 2 RL )i refers to the contributions originated by the analytically contin- 
uation into superspace of the coupling constants of the superpotential, while {m 2 RL )2 to the usual term 
proportional to the [i parameter. Notice that in the sneutrino sector only the LL block is present and that 
any contribution to the sneutrino masses associated to w u can be safetely neglected. 

In order to compute the contributions to the LFV transitions, it is first necessary to move to the physical 
basis, where the kinetic terms are canonical and the charged lepton mass matrix is diagonal. Performing 
the same transformations in all the components of the interested chiral supermultiplets, to avoid flavour- 
violating gaugino-lepton-slepton vertices, it is possible to recover the complete set of non-universal BCs 
at the scale Af. Since the radiative lepton decays occur at lower-energies than A/, it is necessary to evolve 
the soft SUSY breaking terms down to these energies, considering the non-universal BCs. 

We will specify the flavour-dependent contributions for the slepton mass matrices in the physical basis 
in the following, when entering more in the details of the different classes of models, while here we report 
the supersymmetric F— and D— term contributions and (m 2 RL )2' 

{™>I ll )f = mjrrii (m 2 eLL ) D = (~ + sin 2 9w) cos 2/3 m| x 1 (64) 

,2 \ _ n C™2 * ' 1 



«ll)f = (mi LL ) D =\+-\coB2Pm%xl (65) 

{m 2 eRR ) F = m e mj ( m eflii) d = — sin 2 9 W cos 2/3 m| x 1 (66) 

( m e«i)2 = -/xtan/3 m e . (67) 

5.2.2 The LFV Transitions 

Having defined the procedure to get the low-energy expressions for the soft SUSY breaking terms, we can 
now proceed with the radiative lepton decays, \i — > ey, t — > fiy and t — > ej. In this section we provide the 
basic formulae for the computation of the normalized branching ratio Rij for these LFV transitions in the 
so-called mass insertion (MI) approximation ]127[ (for the complete formulae see for example Refs. ]128f 
|130| ). In this approximation, the observables are expressed in terms of the off-diagonal elements of the 
slepton mass matrices msuSY, normalized to their average mass, and the normalized branching ratio 
can be expressed as: 

R ^ = J 8 S a (i^i 2 +i^i 2 )- (68) 

At the LO, the amplitudes A^ and A l R are given by: 

av IX \ i m SUSY rx 
A£ = a LL (d l j) L L + aRL \°ij)RL 

m mi (69) 

Aij fX ^ i m SUSY x 
A R = a RR(°ij)RR + a LR \Oij)LR 
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with ace (C, C — L, R) dimensionless functions of the SUSY parameters. We list in Tab.|4]the expres- 
sions and the numerical values of the functions ace 7 m the limit \i = M\$ = msuSY- The parameters 
See are the mass insertions of the m 2 eCC , : 

[6 ij ) ccl =\^ ) ^. (70) 

" l SUSY 



Table 4 Coefficients dec' characterizing the transition amplitudes for /i — > ey, r — > cy and r — > /try, in 
the MI approximation, taking the limit the /i = Mi, 2 = msuSY and universal BCs. Numerical values are 
given in units of g 2 /(192-k 2 ) and for tan /3 = (2,20). 



aLL 


1 g 2 


1 - 15 tan 2 6 W +4(4 + 


5 tan 2 tan /3 


+(2, 22) 


240 16tt 2 


O-RL = a LR 




12 16^ t&V?9w 




0.30 


iRR 




60lL» tBn9 ^ [ - 6 - 


tan /?] 


-(0.5, 1.6) 



It is interesting now to discuss the main sources of flavour violation that can affect these observables. 
We first consider the case without RH neutrinos, for which the RGE running has a negligible effect to the 
low-energy expressions of the soft SUSY breaking terms, as it is discussed in the App.fA] In this case, it is 
easy to identify the source of flavour violation and to estimate its contributions. 

ll,rr- Contributions to the off-diagonal entries of (Sij)LL arise in correspondence of non-canonical kinetic 
terms due to multiple flavon insertions. Indeed, in the Kahler potential, it is always possible to intro- 
duce a term as with $ a generic flavon, and typically this term produces off -diagonal entries of 
the order of V EV 2 / A 2 . The off -diagonal entries of (S^ll are generically proportional to the same 
quantities. Similar comments apply to the (Sij)mi term, but its impact on the ratios R^ is subdom- 
inant with respect to (Sij)LL, as follows from Tab. |4] In some cases, when the flavor $ transforms 
only under S4 and in particular is a singlet under the additional Abelian symmetries that define G/, it 
is possible to generate d = 3 terms in the Kahler potential. As a result the off-diagonal entries could 
be proportional to VEV/Af. This situation is indeed realized in the models in Refs. 188,89 107 1, 
dealing with the BM pattern. 



rl,lr- The off-diagonal entries of {5ij)RL (and (S.^lr) in the flavour basis have the same origin of the 
charged lepton Yukawa entries. In the mass basis, it is unlucky a change of the order of magnitudes of 
the off-diagonal entries: indeed, this can happen only in the case of precise alignments in the flavour 
space of the contributions to the Yukawa matrix and the slepton mass matrix. The flavour structure of 
the charged lepton Yukawa, and of the RL block of the slepton masses, strongly depends on the class 
of models. Notice that the contributions associated to (Sij)nL are dominant over those associated to 
(Sij)LR, due to a suppressing factor mj/rrii. 

Whether the (Si j)ll contributions or the (Sij)nL ones are the dominant one in the ratios Rij depends on 
the magnitude of VEV/Af and of the off -diagonal entries of the charged lepton Yukawas: both of these 
aspects are model dependent and we will enter in the details of the different class of models in the next 
section. 

When RH neutrinos are present in the spectrum then the direct connection between the contributions 
in the Lagrangian and the flavour violating effects is partially lost. Indeed, from closer look to the RGEs 



(eqs. (88 i-(90i), it follows that large (S^ll are generated even in the case of universal BCs, due to the 



O(l) entries of the Dirac neutrino Yukawa. Depending on the specific lepton mixing and on the flavour 
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BCs of a model, 1) it is possible to neglect the specific flavour BCs and (8ij)i,L is determined only by the 
RGE evolution with universal BCs, or 2) (Sij)LL results from the RGE evolution considering the flavour 
BCs. We will discuss further these possibilities in the next section, when specifying the class of models. 



5.2.3 LFV in TB Models 

For models with the TB mixing, as discussed in Sec. [2j focussing only to the scenario with the highest 
success rate and only on the NH case, when the parameter 9 takes the value of 0.086 then the success rate 
to get the mixing angles inside their present 3<t ranges is maximized. For such value the contributions to 
flavour violating radiative lepton transitions are indeed relevant. 

In order to illustrate the general predictions of the TB models, we assume the following mass matrices 
for the sleptons, in the physical basis: 



( m2 (e,u)LL)K 
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(71) 



(72) 



A 



(73) 



where all the parameters should be consider as random complex numbers with absolute value 0(1), apart 
from no and n| that are positive to favour positive definite square-masses and to avoid electric-charge 
breaking minima and further sources of electroweak symmetry breaking. 

The LL contributions provide the dominant effects on the branching ratio of the lepton radiative transi- 
tions for larger values of the tan (3, while for small values also the RL contributions turn out to be relevant. 
Furthermore, in absence of RH neutrinos, the RGE evolution do not give relevant effects to the low-energy 
expressions of the slepton mass matrix. As a result, entering these expression for the slepton mass matrices 



in eq. (69 1, we find 



Ri 



(\a LL 0(9 2 



(74) 



487r 3 a , , . „ _ „ 

\o-rlO{9)\ 

'F" l SUSY 

for all the three transitions. 

When considering the present bound on BR(u — > erf) from the MEG collaboration, we can constrain 
the supers ymmetric parameter space of the model. The results are shown in Fig. [6] that represents the 
dependence of BR(fi —> e-f) on Mx/2 f° r fixed values of ttt-o- The value for 9 in these plots is 0.103, that 
maximizes the success rate when the corrections to the TB pattern arise from the charged lepton sector, 
which has the best statistical chances. Only the NH case is shown in the plots. The results for the IH case 
as well as those when 9 = 0.104 are very similar. 

When considering models with RH neutrinos, the analysis is slightly different and it is necessary to 



account for the RGE contributions. Considering the flavour BCs in Eqs. (76 1-(78 I, we find that the impact 
on the RG evolution into all the three branching ratios is relevant. In particular, at low-energy, the off- 
diagonal entries of (m 2 LL ) K are multiplied by a factor 3(2), for the NH (IH) case, with respect to the 
corresponding entries at the high-energy. However, this does not change the prediction that 



R 



fi e 



(75) 
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Fig. 6 TB Models. Scatter plots of BR(fi -> erf) as a function of m 1/2 , for tan/3 = (2,20), 9 = 0.103 and 
different values for mo. The horizontal lines show the current MEG bound BR(fi — > ej) < 2.4 x 10 -12 . 



and also the the plots in Fig.[6]can safely apply in the case with RH neutrinos. 



5.2.4 LFV in BM Models 



In this class of models, the lepton mixing matrix at LO coincides with the BM pattern and larger corrections 
are necessary to bring the solar angle in agreement with the data. In the models presented in Refs. [88} 
106 107) , such corrections arise from the charged lepton sector and are of the order of 9 ss O(0. 1), as 



we have previously discussed. No model has been presented so far in which the corrections arise from the 
neutrino sector and therefore we disregard this possibility in the following analysis. 



In this context and considering for the next analysis the specific case of Ref. [ 88 1, the particular charge 
assignment of the flavons and their vacuum alignment determine an hierarchy between the off-diagonal 
entries of the slepton mass matrices. Moreover, a triplet of S4 not transforming under the additional 
Abelian symmetries of Gf is present in the model and as a result it originates non-renormalisable terms in 
the Kahler potential suppressed by only one power of the cut-off scale Af. More in details, in the physical 
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basis we get 
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where all the parameters are treat as in the previous case with TB models. 

The LL contributions provide the dominant effects on the branching ratio of the lepton radiative transi- 
tions. Also for this case, in absence of RH neutrinos, the RGE evolution do not give relevant effects to the 
low-energy expressions of the slepton mass matrix. As a result, entering these expression for the slepton 



mass matrices in eq. ( 69 1, we find 
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As for the previous case with the TB mixing, the most stringent constraint come from the BR(/i —> ej) 
and the results are shown in Fig. IT] The plots represent the dependence of BR(/i — > ej) on M X j 2 for fixed 
values of m . The value for 9 in these plots is 0.156, that maximizes the success rate when the corrections 
to the BM pattern arise from the neutrino sector, which has the best statistical chances. Only the NH case 
is shown in the plots. The results for the IH case as well as those considering 9 = 0.157 are very similar. 

When considering models in which RH neutrinos enrich the spectrum, the analysis can be different. 
Considering the flavour BCs in Eqs. ( 76 H 78 1, we find that the impact on the RG evolution into i? M e and 
R T e is negligible and that only into R T „ is relevant. In particular, at low-energy, the off-diagonal entries of 



L eLL) K 



are all of order of 0(9). As a result the previous prediction for the branching ratios in eq. ( 80 » 



is modified and considering the presence of RH neutrinos it is given by: 



R 



fi e 



(81) 



Considering more in detail the ^ — > ej transition, the plots in Fig. [7] also apply in the case with RH 
neutrinos. 



6 Conclusions 

In the last decade, numerous flavour models based on discrete symmetries have shown their ability in 
accommodating the lepton mixings. In particular, special mixing patterns, that predict the mixing angles 
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M 1/2 (GeV) Mi/2 (GeV) 

(c) m = 200 GeV, tan /3 = 20. (d) m = 5000 GeV, tan (3 = 20. 

Fig. 7 BM Models. Scatter plots of BR(fi -> ej) as a function of m 1/2 , for tan/? = (2,20), 9 = 0.156 and 
different values for mo. The horizontal lines show the current MEG bound BR(fi — > e-y) < 2.4 x 10 -12 . 



independently of the neutrino spectrum, are achieved at the LO: these patterns are the TB, the BM, the GR 
and some TM schemes, that we refer to as FD patterns. 

Not with all the discrete symmetries commonly used it is possible to describe all these special mixing 
patterns. However, the group S4 has been implemented in different models in order to describe both the 
TB, the BM and the TM schemes. In this paper, we have reviewed the main features of the S4 discrete 
group, summarized on the existing models based on this group and analyzed general phenomenological 
signatures of these models. 

A relevant question that we accounted for is about the naturalness of the FD patterns, after that the 
reactor angle has been proven to be non- vanishing. In particular, we analyzed the success rates of TB and 
the BM patterns with general corrections arising either from the charged lepton sector or by the neutrino 
sector. The results favour the TB scheme with general corrections from the charged lepton sector. For this 
case, the success rate reaches about the 15%. Better results could be achieved considering specific type 
of corrections, that for example affect only two of the mixing angles and not all of them. However, in 
this case, the success of the model is due to a dynamical trick, consisting in reducing as much as possible 
the number of NLO independent corrections, that, moreover, must act in specific directions in the flavour 
space. 

Apart the naturalness requirement, a complete model should undergo to a series of criteria: account 
for a mechanism for the flavour symmetry breaking; describe the quark sector in a unified approach, i.e. 
using the same flavour symmetry both for quarks and leptons; produce predictions in the lepton and quark 
sectors. More criteria are satisfied more a model should be considered complete even if it is true that nature 
could have been chosen an unexpected way. 
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In order to test these models and to distinguish one model from the others, it is necessary to investigate 
on their phenomenological signatures. We pointed out that often sum rules are present in discrete flavour 
models and as a result quite clear correlations arise among the sum of the neutrino masses, the kinematic 
electron mass in the beta decay and the neutrino-less double beta decay effective mass. Alternatively, lepton 
flavour violating transitions, such as £^ — > £,-7, provide strong constraints, once assuming the presence of 
flavour-sensitive new physics at the TeV scale. In this paper, we concentrated in the supersymmetric 
case, distinguishing the case in which the neutrino masses are explained through the Weinberg operator or 
through the type I See-Saw mechanism. Even if the contributions from the RGE are different in the two 
cases, the constraints arising form the /i — > ej decay are similar: these models are not compatible with 
very light supersymmetry. 

Even though the huge effort of these years in constructing flavour models to describe masses and mix- 
ings for the neutrinos, and more in general for all the fermions, it is discouraging that no illuminating 
strategy arise form this scenario. On the other hand, this is partially related to the large uncertainties still 
present in the flavour sector. The hope is that with a better determination of the lepton mixing angles and 
with the knowledge of the CP phases, the neutrino mass scale, the type of the neutrino nature and spectrum, 
it will be finally possible to shed light on the origin of the fermion masses and mixings. 
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- Institute for Advanced Study, funded by the German Excellence Initiative. 



A Renormalisation Group Effects 



The renormalization group equations in the MSSM context for the soft mass terms (m 
and A e = y/2( 

m "eRL)i/ iv cos @)> denoting t' = log(A^,/ tiisusy), are [ |128||l3l| : 
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(84) 



1^,-Y eij - 

(85) 

where gi^ 2 are the gauge couplings |^] of SU(2)l x U(1)y, Mi,2 the corresponding gaugino mass terms, 
Y e ,d = tnt,d I v d are the Yukawa matrices for charged leptons and down quarks, Ad — (tn 2 RL )i / Vd and: 



S = Tr (m 2 qLL 



dRR 



2m lRR 



{ml LL ) K + (m 2 eRR ) K ) - m 2 Hd + m 2 Hu 



(86) 



In the GUT normalization, such that Q2 = g and g\ 
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The matrix (m 2 LL )K coincides with {m 2 eLL )K and has the same evolution. For squarks we have intro- 
duced soft mass terms analogous to those previously discussed for sleptons. To estimate the corrections to 
the slepton masses induced by the renormalization group evolution we adopt the leading logarithmic ap- 
proximation and substitute each of the running quantities with their initial conditions at the scale ~ A/ 



in eqs. (82i-(85i 



In this approximation and considering that in general the off-diagonal entries of the soft mass matrices 
are smaller than the corresponding elements into the diagonal, one easily sees that the largest corrections to 
the matrices (m 2 e u \ ll )k an d ( m eRR)K come from electroweak gauge interactions and are proportional 
to the identity matrix in flavour space. Due to the negative sign of the dominant contribution these diagonal 
elements increase by evolving the mass matrices from the cutoff scale down to the electroweak scale. 

The largest corrections on the off-diagonal elements of the LL and RR low-energy soft mass are in the 
(13) and (23) sectors: 



{ m2 (eM)LL)K i3 
( m lRR)K i3 



7^2 Vr (m (eiV)LL )K i3 log I 



1 



I67P 



Vr {™l RR )K i3 log 



msusY 

—) 

™ SUSY J 



(87) 



All such contributions can be safely neglected, because they are suppressed by y T and by log (Al/itisusy)- 
The matrix A e gets a first correction by an overall multiplicative factor that can be absorbed, for in- 
stance, by a common rescaling of the parameters, and a second correction of the type A e — > A e + S Y e , 
which is only into the diagonal elements and easily reabsorbable redefining the parameters. Finally, such 
corrections to the off-diagonal elements of A e are negligible. We can conclude that the corrections on the 
off -diagonal entries of the soft mass matrices induced by the RG running are either negligible or could be 
absorbed into the parametrization. 

When considering the See-Saw case, the RGEs that describe the running of the sfermion masses are 
different | fT28l[T3Tl : 



d/i 



d/i 
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™ZlL 
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16tt 2 



(88) 

(89) 
(90) 



where Y v is the neutrino Dirac Yukawa in the basis of diagonal charged leptons and RH neutrinos. Since 
two of the lepton mixing angles are large, then Y v in general has large off-diagonal entries. This would 
contribute to flavour violating processes much more than in the case without RH neutrinos, especially from 
the LL sector. 

The low-energy quantities should be calculated considering the flavour boundary conditions at the scale 
A^ ss Af and evolving them down by the use of eqs. (88 l-(90l. This can easily be done numerically, but 
not analytically. 
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